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Abstract 

Ages of the magnetar IE 2259+586 and the assoeiated supernova remnant CTB 109 
were studied. Analyzing the Suzaku data of CTB 109, its age was estimated to be ~14 kyr, 
which is much shorter than the measured characteristic age of IE 2259+586, 230 kyr. This 
reconfirms the previously reported age discrepancy of this magnetar/remnant association, 
and suggests that the characteristic ages of magnetars are generally over-estimated as com¬ 
pared to their true ages. This discrepancy is thought to arise because the former are cal¬ 
culated without considering decay of the magnetic fields. This novel view is supported 
independently by much stronger Galactic-plane concentration of magnetars than other pul¬ 
sars. The process of magnetic field decay in magnetars is mathematically modeled. It 
is implied that magnetars are much younger objects than previously considered, and can 
dominate new-born neutron stars. 
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1. INTRODUCTION 


1.1. Relations between Magnetars and SNRs 

Twenty seven Galaetie and Magellanie X-ray sources are thought to constitute a class of ob¬ 
jects called magnetars, which are single neutron stars (NSs) with extremely strong magnetic fields of 
B = G. They are believed to shine (mainly in X-rays) consuming the energies in their strong 

magnetic fields, because their X-ray luminosities exceed their spin-down luminosities and they are 
not likely to be accreting objects. The magnetar concept well explains other peculiar characteristics 
of these objects, such as long pulse periods clustered in a narrow range (2 — 11 s), relatively large 
period derivatives, and unpredictable sporadic burst activities. However, we do not know yet how 
they are formed and how such strong fields evolved. 

Supernova remnants (SNRs) associated with magnetars are expected to provide us with valu¬ 
able clues to the scenario of magnetar production (e.g., Vink 2008; Safi-Harb & Kumar 2013). As 
a result, the study of SNR/NS associations (e.g., Seward 1985; Chevalier 2005, 2011) has been re¬ 
activated since the concept of magnetar has emerged. Although no clear difference in the explosion 
energy has yet been found between SNRs with and without magnetars (e.g., Vink & Kuiper 2006), 
an X-ray metallicity study of the SNR Kes 73, hosting the magnetar IE 1841-045, led Kumar et al. 
(2014) to infer that the progenitor of this system had a mass of ^ 20 Mq, where Mq is the solar mass. 
Through investigations of several SNR/magnetar associations, Safi-Harb & Kumar (2013) character¬ 
ized environments that are responsible for the magnetar production, and reinforced the view of rather 
massive progenitors. 

Apart from the progenitor issue, one particularly interesting aspect of magnetars, which can 
be studied by simultaneously considering the associated SNRs, is their age comparison. Of course, 
as discussed by Gaensler (2004), the ages of a magnetar and of the associated SNR, estimated in¬ 
dependently, must agree for them to be regarded as a true association. However, these two age es¬ 
timations sometimes disagree even in pairs with very good positional coincidence, including the IE 
2259-1-586/CTB 109 pair which is the topic of the present paper. This is often called “age problem”. 
While Allen and Horvath (2004) suspected that the problem arises because an SNR age estimate is 
affected by the presence of a magnetar, Colpi et al. (2000) instead attributed it to magnetic field decay 
of magnetars, which can make their characteristic ages much longer than their true ages. After the 
discovery of SGR 0418-1-5729 (Rea et al., 2010), a magnetar with a low dipole magnetic field, the 
field’s decay scenario has become more attractive (Dall’Osso et al., 2012; Igoshev, 2012). 

Through X-ray studies of CTB 109, the present paper attempts to address two issues. One is to 
solve the age problem with the IE 2259-1-586/CTB 109 system, and the other is to conversely utilize 
the result to reinforce the nature of magnetars as magnetically driven NSs. After a brief introduction 
to the target system (section 2), we describe in section 3 and section 4 recent Suzaku observations 
of CTB 109, and reconfirm the age problem in the system (section 5). Then, an attempt is made in 
section 6 to solve it invoking the decay of magnetic fields. Einally, we discuss some implications 
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for the general view of NSs, ineluding in partieular their magnetie evolution. Other topies will be 
discussed elsewhere, including more detailed X-ray diagnostics of CTB 109, the origin of its peculiar 
half-moon shape, and characterization of the progenitor. 

2. Magnetar IE 2259+586 and SNR CTB 109 

2.1. CTB 109 

The Galactic SNR, CTB 109, hosting the central point X-ray source IE 2259+586, was first 
discovered by the Einstein Observatory (Gregory & Eahlman, 1980) as an extended X-ray source 
with a peculiar semi-circular shape. It was independently identified as a shell-type SNR by radio 
observations at 610 MHz (Hughes et ah, 1981). A 10 GHz radio map taken with the Nobeyama Radio 
Observatory revealed good positional coincidence between the radio and X-ray shells, while detected 
no significant enhancement from IE 2259+586 (Sofue et ah, 1983). 

Through CO molecular line observations, Heydari-Malayeri et al. (1981), Tatematsu et al. 
(1985), and Tatematsu et al. (1987) found a giant molecular cloud located next to CTB 109, and 
suggested that it may have disturbed the SNR on the west side. Sasaki et al. (2004) conducted a 
comprehensive X-ray study of this SNR with XMM-Newton. Assuming a distance of 79 = 3.0 kpc 
(Kothes et al., 2002; Kothes & Eoster, 2012), they estimated the shock velocity, age, and the explosion 
energy as Vs = 720 ± 60 km s“\ 8.8 kyr, and (7.4 ± 2.9) x 10^° erg respectively. Strong evidence for 
an interaction between the SNR shock front and the CO cloud was found by using ^^CO, ^^CO and 
Chandra observations (Sasaki et al., 2006). Eurthermore, gamma-ray emission was detected with the 
Eermi-Lat from CTB 109 (Castro et al., 2012). Einally, using Chandra, Sasaki et al. (2013) detected 
emission from the ejecta component and refined the age as 14 kyr. 

2.2. IE 2259+586 

The compact object IE 2259+586 was first detected in X-rays nearly at the center of CTB 109 
(Gregory & Eahlman, 1980). It was soon found to be a pulsar, and the pulse period was at first 
considered as P = 3.49 s (Eahlman et al., 1982). This was due to the double-peaked pulse profile, and 
the fundamental period was soon revised to P = 6.98 s (Eahlman & Gregory, 1983). Repeated X-ray 
observations enabled the spin down rate to be measured as P = (3 — 6) x 10“^^ ss“^ (Koyama et al., 
1987; Hanson et al., 1988; Iwasawa et al., 1992), and these results made it clear that the spin-down 
luminosity of IE 2259+586 (5.6 x 10^^ erg s“^) is far insufficient to explain its X-ray luminosity, 
1.7 X 10^"^ erg s“k Due to this and the long pulse period, IE 2259+586 was long thought to be an 
X-ray binary with an orbital period of ~ 2300 s (e.g. Eahlman et al. 1982), and extensive search for 
a counterpart continued (Davies et al., 1989; Coe & Jones, 1992; Coe et al., 1994). However, no 
counterpart was found (Hulleman et al., 2000), and instead, tight upper limits on the orbital Doppler- 
modulation have been obtained as axSini < 0.8 light-s (Koyama et al., 1989), a^sini < 0.6 light-s 
(Mereghetti et al., 1998), and a 3 ;Sini < 0.028 light-s (Baykal et al., 1998). Here, is the semi-major 
axis of the pulsar’s orbit, and i is the orbital inclination. These strange properties of IE 2259+586 led 
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this and a few other similar objeets to be ealled Anomalous X-ray Pulsars (AXPs). 

In the 1990’s, several attempts were made to explain IE 2259-1-586 without invoking a eom- 
panion: e.g., massive white dwarf model (Usov, 1994), and preeessing white dwarf model (Pandey, 
1996). Monthly observations of IE 2259-1-596 over 2.6 years with RXTE gave phase-eoherent timing 
solutions indieating a strong stability over that period (Kaspi et ah, 1999). This favored non-aeeretion 
interpretation. Heyl & Hemquist (1999) suggested that spin-down irregularities of AXPs are sta- 
tistieally similar to glitehes of radio pulsars. Meanwhile, the eoneept of magnetar was proposed to 
explain Soft Gamma Repeaters (SGRs) as magnetieally-powered NSs, namely, magnetars (Dunean 
& Thompson, 1992; Thompson & Dunean, 1995). Eurthermore, like SGRs, IE 2259-1-586 showed 
an X-ray outburst (Kaspi et ah, 2002; Gavriil et ah, 2004; Woods et ah, 2004). Today, AXPs inelud¬ 
ing IE 2259-1-586, as well as SGRs, are both eonsidered to be magnetars. Employing the eanonieal 
assumption of spin down due to magnetie dipole radiation, the measured pulse period of P = 6.98 s 
(Iwasawa et ah, 1992), and its derivative, P = 4.8 x 10“^^ss“^ (e.g., Gavriil & Kaspi 2002), give a 
dipole magnetie field of 5.9 x 10^^ G and a eharaeteristie age of Tc = P/2P =230 kyr. 

3. OBSERVATIONS AND DATA REDUCTION 

The Suzaku observations of IE 2259-1-586 and CTB 109 were made on two occasions. The 
first of them was conducted on 2009 May 25 as a part of the A04 Key Project on magnetars (Enoto 
et ah, 2010). The three cameras (XISO, XISl, and XIS3) of the X-ray Imaging Spectrometer (XIS) 
onboard Suzaku were operated in 1/4-window mode with a time resolution of 2 s (Koyama et ah, 
2007), to study the 6.98 s pulsation. The rectangular (17' x 4'.3) fields of view of the 1/4-window 
were center on this magnetar, while the SNR was partially covered. As the second observation, we 
performed new four pointings onto CTB 109 with Suzaku (PLT.Nakano) on 2011 December 13 . In 
order to synthesize a whole image of the SNR, all XIS cameras were operated in full-window mode 
under the sacrifice of time resolution. A log of these observations is given in Table 1. 

In the present work, we use the XIS data which were prepared with version 2.7.16.31 pipeline 
proceeding, and the calibration data updated in January 2013. The data of the Hard X-ray Detector 
(HXD) from neither observation are utilized here, since the Suzaku data analysis in the present paper 
focuses on CTB 109, rather than the central magnetar. The data reduction was carried out using the 
HEADAS software package version 6.13, and spectral fitting was performed with xspec version 
12.8.0. The redistribution matrix files and the auxiliary response files of the XIS were generated with 
xisrmfgen and xissimarfgen (Ishisaki et ah, 2007) respectively. 

4. DATA ANALYSIS AND RESULT 

4.1. Image Analysis 

Eigure 1 shows a gray-scale mosaic image of CTB 109 obtained with the Suzaku XIS, shown 
after subtracting non X-ray background. The magnetar, IE 2259-1-586, appears as a bright source at 


Table 1. Log of Suzaku observations of CTB 109. 


Observation ID 

a 

<5 

Start Time 

Exposure) ks) 

404076010 

23‘'01“04".08 

58°58T5".6 

2009-05-25 20:00:17 

122.6 

506037010 

23‘'01“06".96 

59°00'50".4 

2011-12-13 06:48:41 

40.8 

506038010 

23^ 00“ 26k88 

58°44T3".2 

2011-12-14 04:47:02 

41.4 

506039010 

23^ 03“ 06k96 

58°58'51".6 

2011-12-15 01:57:25 

30.4 

506040010 

23^ 03“ 06k96 

58°40'51".6 

2011-12-15 18:03:52 

30.5 


the center. The crisscross region including IE 2259+586 was taken in the first observation, while 
four square regions represent those from the second one. Thus, the mosaic XIS image reconfirms the 
half-moon like morphology of this SNR (section2.1). 

As mentioned in section 2.1, the lack of western part of CTB 109 is usually attributed to 
its interactions with giant molecular clouds. Bright regions in the SNR are also thought to be the 
signature of such an interaction (Sasaki et ah, 2006). These issues, together with detailed spatial 
distributions of X-ray spectral properties, will be postponed to another publication. The present paper 
instead deals with average X-ray spectra, because our prime motivation is to estimate the age of 
CTB 109. 

4.2. Spectral Analysis 

Figure 2 shows XISO, XISl and XIS3 spectra of CTB 109 from the second observation, ex¬ 
tracted from the regions indicated in figure 1 as “NE” and “SE”, which represent two eastern pointing 
positions. Since these eastern parts of the SNR has kept smooth round shape, the effect from the inter¬ 
action with CMC seems to be small. Background spectra were extracted from source free regions in 
the same observation which are also indicated as “BCD” in figure 1. Spectral bins were summed up 
to attain a minimum of 30 counts bin“^. The ”NE” and ”SE” spectra are very similar and both exhibit 
emission lines due to highly ionized atoms such as Ne-IX triplet (~0.92 keV), Ne-X Eya (1.02 keV), 
Mg-XI triplet (~1.35 keV), Si-XIII triplet (~1.87 keV), and S-XV triplet (~2.45 keV). 

We first applied a variable-abundance non-equilibrium ionization (VNEI) plasma emission 
model to the NE spectra in figure 2 (a). However, even when abundances of Ne, Mg, Si and S 
are allowed to vary freely, the reduced of fho At was not made lower than 1.5. Thus, the single 
temperature VNEI model was rejected. Other plasma models in xspec, such as vpshock and 
vmekal were also unsuccessful. 

Then, we considered that the SNR emission consists of two components including ejecta and 
shocked inter-stellar medium (ISM), and added a non-equilibrium ionization (NEI) plasma emission 
model to model shocked ISM component. Abundances of ISM component was fixed to solar to reduce 
the number of free parameters. 

By introducing this two-component emission model, the fit was improved to jv =1.28. 
Since different XIS cameras gave discrepant fit residuals around Mg-XI Ka and Si-XIII Ka lines. 
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Fig. 1. A mosaic image of CTB 109 obtained in 0.4-5.0 keV with the Suzaku XIS. After sub¬ 
tracting the non X-ray background, the image was corrected for exposure and vignetting. The 
on-source and background regions are indicated by white and dashed black lines, respectively. 
Some corners of the square XIS fields of view are masked to remove calibration isotopes. 

presumably due to ealibration uneertainties of the XIS, we allowed gain parameters of the XIS eam- 
eras to vary. Then, the fit has beeome acceptable with jv = 1.06 (1293/1223). The obtained best-fit 
parameters are shown in table 2, while the gain parameters in table 3. Since the obtained abundances 
of the second (Plasma 2) component are all consistent with 1 solar, it is also likely to be dominated by 
the ISM. The two components are both inferred to be somewhat deviated from ionization equilibrium. 
The fit was not improved even when the abundance of the Plasma 2 component is allowed to change 
from 1.0. Similarly, we analyzed the SE spectra shown in figure 2 (b), and derived the parameters 
which are again summarized in table 2. Thus, the two temperatures are in good agreement between the 
two regions, although the SE region gives somewhat higher metal abundances. This two-component 
emission model was also employed by Sasaki et al. (2013). 

5. Estimation of the Age of CTB 109 

Now that the average plasma properties of CTB 109 have been quantified, let us proceed to 
our prime goal of studying this SNR, i.e., its age estimation. Eirst, the physical radius of CTB 109 
is estimated as R = (16 ± l)d 3,2 pc , from its angular size of ~ 16' and the estimated distance 
D = 3.2 ± 0.2 kpc (Kothes & Eoster, 2012). Here, ^ 3.2 = D/3.2 is scale factor of distance. Next, 
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Table 2. Best-fit spectrum pai'ameters for the NE and SE regions*. 


Component 

Parameter 

NE 

SE 

Absorption 

Ah (lO^^ cm“^) 

0.78 ±0.01 

n 7Q+0.01 
' *^-0.02 

Plasma 1 (VNEI) 

kTi 


0 (^15+0.02 

'^•'^^-0.01 


rji* (lO^^ cm“^ s) 

n 07+0.05 
-0.04 

n OO+0.05 
'^■^^-0.06 


Ne (solar) 

0.8T0.2 

0.8+? 


Mg (solar) 


I.ITO.I 


Si (solar) 

1.2±0.1 

1 Y+0.1 

^■‘-0.3 


S (solar) 

l.OTO.l 

1 7+0.2 
-^■'-0.3 


Fe (solar) 

Q QQ+0 05 
^■^^-0.06 

-*-■^-0.1 


Alt (10-2 cm-5) 

Q 7+0.2 
-0.1 

1.2±0.1 

Plasma 2 (NEI) 

m 

0.26 ±0.01 

0.25 ±0.01 


T] 2 * (10^2 cm-3 s) 

> 1 

> 1 


Abundance (solar) 

1 (fixed) 

1 (fixed) 


Aat (10-2 cm-5) 

19.6+;^ 

13.2+;® 

X^/dof 

wn -:r“!— ir-r: -rr 


1294/1223 

1138/951 


* Uncertainties are statistical errors at 90 % confidence. 

* Ionization parameter, defined as p = rief 


t Normalization of NEI or VNEI, defined as K = f ngnudV 


Table 3. The best-fit gain parameters*. 


region 

instrument 

slope 

offset (eV) 

NE 

XISO 

0.997 ±0.001 

-5.3±0.2 


XISl 

1.004 ±0.001 

-3.0 ±0.3 


XIS2 

0.986 ±0.001 

-7.0 ±0.4 

SE 

XISO 

1.003 ±0.001 

-10.6±0.1 


XISl 

0.993 ±0.001 

6.6±0.1 


XIS2 

0.985 ±0.001 

6.6±0.2 


* Uncertainties are statistical errors at 90 % confidence. 


CTB 109 may be considered in the Sedov phase (neglecting the missing western part). Then, applying 
the Sedov-Taylor similarity solution (Sedov, 1959; Taylor, 1950) to this SNR, its age can be obtained 
as 

2R 

"rSNE = 7— (1) 

where Vg again represents the shock front velocity. 

Assuming the strong shock limit, we can calculate Vg from the post-shock temperature Tpg as 
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Fig. 2. Background-subtracted XISO (black), XISl (red) and XIS3 (green) spectra of CTB 109, extracted from the 
NE (panel a) and SE (panel b) regions. They are fitted simultaneously with a two-component model described in 
the text. The dotted lines indicate individual model components, with the three colors specifying the three cameras. 


V. = \ ^kT, 


3m 


ps 


( 2 ) 


where k and fh, respeetively represent the Boltzmann eonstant and the mean mass per free partiele. 
Assuming a solar abundanee (seetion 4.2), we employed m ~ 0.61 mp where mp is the proton mass. 
Since CTB 109 is a middle-aged SNR without too strong shocks, the electron temperature of the NEI 
component measured in section 4 can be considered to be close to the kinematic ISM temperature, and 
hence to Tps (Ghavamian et ah, 2007). Thus, substituting 0.25 ± 0.02 keV for Tps, equation (2) gives 
Vs = 460 ±40 km s“^, and then equation (1) yields tsnr — (14 ± 2)ds,2 kyr in agreement with Sasaki 
et al. (2013). Compared to the estimated tsnr, the characteristic age of IE 2259-1-586, Tc = 230 kyr 
(section 2.2), is ~16 times larger. 

If we assume that CTB 109 is in cooling phase rather than Sedov phase, the time dependence 
of the radius becomes R oc (McKee & Ostriker, 1977). Then, the age estimation slightly changes 
to 


Tsnr = ^ —, (3) 

7 Vs 

and Tsnr = (10 ± l)d 3,2 kyr is obtained. Therefore the age discrepancy still persists (even increases). 

Eet us cross-check the above estimates using the parameter r] = nj, (table 2) of the ejecta 
component. We assume that the SNR has a half spherical shell with thickness of Ai? = R/12 (as¬ 
suming 4/37ri7^?7,omp = AnR?^RuQmp ), and the ISM component corresponding to Plasma 2 in 
table 2 is emitted from this shell. Eurthermore, as a simplest approximation, the ejecta component 
(Plasma 1 in table 2) may be assumed to uniformly fill the inner region, considering that the re¬ 
verse shock has reached the center of the SNR. In the NE region, the extracted emission volumes 
of Plasma 1 (ejecta) and Plasma 2(ISM) can be obtained numerically as Ei = 4.0 x 10^®(i|2 
and V 2 = 1.4 x 10^®(i|2 respectively. Then , the spectrum normalizations in table 2 give 
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the averaged density of the ejecta as rii = (0.33 ± cm“^ and that of the ISM-shell as 

77.2 = (1.2 ± Q.l)d'^y‘^ cm“^. The time required for the ejecta to become ionized as we now ob¬ 
serve is hence estimated as 77 / 77.1 = (29 — 3A)d]J^ kyr. Applying the same argument to the SE region 
having the emission volumes of Vi = 6.7 x 10 ^^(i 3 2 cm^ and V 2 = 3.3 x 10®^d| 2 we obtain 
77.1 = (0.46 ± O.OSjdJ/g'^^ cm“^ and 77.2 = (0.55 ± O.OOjdJ/^'^^ cm“^, and 77/771 = (11 — lOjcig^^ kyr. 
Pre-shock density is estimated as Uq = 772/4 = (0.1 — 0 . 3 )(i^ 2 ^^ cm“^. Even though these estimates of 
t'snr niust have a certain range of systematic uncertainties, large discrepancies significantly remain 
between tsnr and Tc- In fact, it would be difficult to think that CTB 109 are emitting X-rays even 
at an age of 230 kyr while keeping the regular shape (except the missing western half), because its 
density environment as estimated is quite typical of a solar neighborhood. Thus, we reconfirm the 
previously reported age discrepancy (Sasaki et ah, 2013) between CTB 109 and IE 2259-1-586. 

6. Solving the Age Discrepancy 


The magnetar IE 2259-1-586 is located nearly at the very center of the half-moon-shaped shell 
of CTB 109 (section 4.1, figure 1). This coincidence is difficult to explain by invoking a chance 
superposition of the two objects. We therefore assume that IE 2259-1-586 and CTB 109 were in¬ 
deed produced by the same supernova explosion (section 1.1), while Tc of IE 2259-1-586 is somehow 
significantly overestimated, compared to its true age which we consider to be close to tsnr- 

6.1. Case with a Constant Magnetic Field 


To solve the issue of the suggested overestimate of Tc after Colpi et al. (2000) and Dall’Osso 
et al. (2012), let us begin with reviewing the meaning of t^. In general, the spin evolution of a pulsar 
with dipole surface magnetic field B is expressed empirically as 

^ (4) 

dt 

with b = 327r^i?pgj./3//ioC^ and the braking index of 77. = 3, where Rp^r = 10 km is the pulsar’s radius, 
/ = 9.5 X 10"^“^ g cm^ its momentum of inertia, /iq vacuum permeability and c the light velocity. If we 
use the pulse period P = 27:/uj and its derivative P instead of uj and uj, equation (4) becomes 


B = 


I pp 


~3.2 X IO^VPPG. 


If B does not depend on time t, equation (4) can be integrated as 


t = 


1 


u 




= Tr. 


1 


(5) 


( 6 ) 


77. — 1 \ 6 j , 

where oj and 6j both refer to the present values, while ojq is the angular frequency at f = 0 (i.e. the 

can be neglected, the characteristic age is defined as 

UJ P 


birth). Assuming that {uj/ujqY ^ 


Tr. = 


(77-1)0; ( 77 ,-l)P 


(7) 
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These equations are generally used for pulsars, and found in some textbooks (e.g., Lyne & Graham- 
Smith 1998). 

More generally, the true age of the pulsar, denoted by to, ean be compared with its Tc as 


I2. 

to 




1 + 



( 8 ) 


where Pq = 27r/a;o, and the last expression is the first-order approximation in (Pq/P)" Thus, 


becomes somewhat larger than to if (Pq/P)"^ cannot be neglected. Conversely, if we somehow have 
an independent estimate of to its comparison with Tc can be used to infer Pq as 


Po = P(-- + l 




(9) 


V Tc y 

For example, the Crab pulsar (Staelin & Reifenstein, 1968), with P = 33 ms, P = 2.42 x 10“^^ ss“^ 
and n = 2.509 (Lyne et ah, 1993), has Tc = 1241 yr. Comparing this with its true age of 960 yr (as 
of 2014), equation (9) yields Pq = 18 ms if assuming n = 2.509, on Pq = 15.7 ms if n = 3.0 (for 
ideal magnetic dipole radiation). Thus, regardless of the employed value of n, the small difference 
between Tc and to of the Crab pulsar can be understood to imply that it has so far lost ~ 3/4 of its 
initial rotational energy in rs-/ 1 kyr. 

In contrast to the above case of young active pulsars, we would need to invoke Pq = P x 0.97 = 
6.76 s, if equation (9) with n = 3 were used to explain the large discrepancy, Tc/to ~ Tc/tsnr ~ 16, 
found in the CTB/1E2259-I-586 system. This would lead to a view that IE 2259-1-586 was born some 
14 kyr ago as a slow rotator of which the spin period is much longer than those (0.2 s to 2 s) of the 
majority of currently observed (hence relatively old) radio pulsars, and has so far lost only a tiny 
fraction of its rotational energy in 14 kyr. However, such a view is opposite to a general consensus 
that new-born magnetars must be rotating rapidly, even faster than ordinary pulsars, in order for 
them to acquire the strong magnetic fields (e.g., Usov 1992, Duncan & Thompson 1996, Lyons et 
al. 2010). Furthermore, an NS with Pq = 6.67 s, has an angular momentum of only 10 ^ of those 
of typical new-born pulsars with Pq ~ 10 ms including the Crab pulsar, and hence would require an 
extreme fine tuning in the progenitor-to-NS angular momentum transfer during the explosion. We 
therefore conclude that the age problem of 1E2259-I-586 cannot be solved as long as its magnetic field 
is assumed to have been constant since its birth. 


6.2. Ejfects of Magnetic Field Decay 

Since the age problem of IE 2259-1-586/CTB 109 cannot be solved as long as B is considered 
constant, we may next examine the case where B decays with time (section 1.1). In fact, the X-ray 
emission of magnetar is thought to arise when their magnetic energies are consumed (Thompson & 
Duncan, 1995). Then, the calculations presented in section 5.1 would be no longer valid, and we need 
to integrate equation (6) considering the time evolution of B. 

Let us consider a simple magnetic filed decay model employed by Colpi et al. (2000), namely 
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dB 

dt 


where a > 0 is a parameter, and a is another positive eonstant. This equation ean be solved as 

Bo 


B(i) = 


(l + f) 

Bq exp 


at) 


Td 


(a ^ 0) 


(a = 0) 


( 11 ) 


where Bq represents the initial value of B, and Td = {l/aB^), an arbitrary constant, means a typical 
lead time till the power-law like decay of B begins. 

Substituting equation (11) into equation (4), we can derive P as a function of t. Then, as 
already given by Dall’Osso et al. (2012), can be expressed as a function of t, Pq , a and Td as 


2 — a 



Tc = 


Td 

2 



Pd/ 


exp 



-1 


(a = 0) (12) 



ro + -ln 



(a = 2) 


where Tq = Po/2Po is the initial value of Tc. The first form of equation (12) reduces to equation (8) 
for a —)■ oo or Td —)■ oo, i.e., the case of a constant B. 


6.3. Magnetic Field Evolution of IE 2259+586 

Our next task is to examine whether the observed values of P and P of IE 2259-1-586 can 
be explained with the picture presented in section 6.2. Equation (12) involves four free parameters, 
namely a, Pq > P and Pq, whereas we have only two observables, P and P at f = to —14 kyr. In section 
5.1, we showed that the effect of Pq can be neglected, when the current P is long enough. Therefore, 
we chose to fix Pq at 3 ms where strong dynamo works efficiently (e.g., Duncan & Thompson 1996). 
To visualize effects of Pq, another solution with Pq = 10,100 ms and a = 1.2 is also shown in figure 
3 [panels (b), (c), and (d)]. Thus, the effects of Pq are limited to very early (<C 1 s) stages of the 
evolution, and its value does not affect our discussion as long as it is much shorter than ^ 6.7 s. 

Then, if a is specified, we can find a pair (Pq^Ti) that can simultaneously explain P and 
P at present. Eigure 3 shows the behavior of such a family of solutions to equation (12). Below, 
we try to constrain the values of a (hence of Bq and tq), assuming that a is relatively common 
among magnetars. This is because the broad-band X-ray spectra of magnetars are determined rather 
uniquely by Tc (Enoto et al., 2010), so that Tc is considered to be tightly related to to even if these two 
are unlikely to be identical: object-to-object scatter in a would cause a scatter in the Tc/Iq ratio, and 
would make the relation of Enoto et al. (2010) difficult to interpret. 
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When a is small (0 < a < 0.5), the field would decay, as seen in equation (10), either exponen¬ 
tially (if a = 0 ) on a time scale of ra, or if a 7 ^ 0 , with a relatively steep power-law after a long lead 
time Td ~ to- The required initial field, Bq ~ 10^® G, is reasonable. However, the implied view would 
be rather ad-hoc: IE 2259-1-586 had been relatively inactive until recently, when it suddenly started to 
release its magnetic energy with a high rate. Furthermore, if such a small value of a were common 
to magnetars, their age differences would make, as in figure 3 (a), their present-day field distribution 
scatter much more widely than is observed. Hence we regard these small values of a unlikely. 

As a increases towards 2.0, the power-law field decay becomes milder, with shorter values 
of Td and stronger initial fields Bq. The implication is that the object started releasing its magnetic 
energy rather soon after the birth, and had already dumped away a large fraction of its rotational 
energy at a very early stage when the field was still very strong. As seen in figure 3 (c), the spin 
period has almost converged to its terminal value (see also Dall’Osso et al. 2012). Therefore, this 
case can explain the observed narrow scatter in P of magnetars, assuming that they share relatively 
similar values of a and Bq. However, the cases with a ~ 2.0 or larger would require too strong initial 
fields, e.g., Bq ~ 10^^ G, which would be much higher than the strongest dipole field observed from 
magnetars, B = 2.4 x 10^^ G of SGR 1806—20 (Nakagawa et al. 2009). Therefore, such large-a 
solutions are unlikely, too. 

To summarize these examinations, figure 3 (e) shows the locus of the allowed solutions on the 
(a,rd) plane, where the values of Bq are also indicated. We thus reconfirm the above considerations, 
that the range of 1 ^ a < 2 is appropriate, in agreement with the suggestion by Dall’Osso et al. (2012). 
Some discussions follow in subsection 7.1. 

7. Discussion 

7.1. Comparison with Other Objects 

We reconfirmed the age problem of IE 2259-1-586 and CTB 109, and presented a way to solve 
it with a simple magnetic field decay model. The result agrees with the basic concept of magnetar hy¬ 
pothesis which implies that the energies stored by their magnetic fields should be consumed to supply 
their X-ray luminosities exceeding those available with their spin down. The amount of released field 
energies can be reflected in the overestimations of the characteristic ages. 

Eet us then examine whether this concept applies to other NS/SNR associations, including 
both magnetars and ordinary pulsars. Figure 4 shows relations between Tc of such single pulsars 
and the ages of their host SNRs. Data points of ordinary pulsars are distributed around the line 
representing Tc/tsnr = 1 with a few exceptions (e.g., Torii et al. 1999 for J1811 —1925/G11.2—0.3). 
Therefore, radio pulsars, including the particular case of the Crab pulsar (section 6.1), are considered 
to be free from the age problem. 

In addition to the ordinary pulsars, figure 4 shows a few other magnetar/SNR associations. The 
magnetar CXOU J171405.7-381031 has a very small characteristic age of 0.96 kyr (Sato et al., 2010), 
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time (year) 



(f) a=0.0,r^=6.6e+03yr,5o=4.9e+14G - 

«=0.6, rj=2.5e+03 yr, 5fl=6.5e+14 G - 

«=1.0,r^=9.2e+02yr,So=9.4e+15G - 

Pg= 10 ms - 

Pq= 100 ms . 

«=1.4,r^=1.6e+02yr,So=1.8e+15G - 

«=2.0, r^=2.5e-03 yr, So=2.0e+17 G - 


Fig. 3. Possible evolution tracks of IE 2259+586 assuming equation (4) and equation (10). Panels (a)-(d) 
represent the behavior of the magnetic field B, the over-estimation factor of the characteristic age (i.e., 
Tc/to), the pulse period P, and its time derivative P, respectively. The six representative tracks are 
all constrained to reproduce the presently measured P and P at t = 14 kyr. The dashed and dotted 
ones assume Pq = 10 ms and 100 ms, respectively, while the other five all Pq = 3 ms. Panel (e) 
shows the trajectory of solutions that can explain the present-day (t = 14 kyr) IE 2259+586. Dashed 
lines indicate the initial filed value Bq. Panel (f) summarize the parameter sets of the trajectories. 

which is consistent, within rather large errors, with the age (O.GSIqJ kyr; Nakamura et al. 2009) the 
associated SNR, CTB 37B. Another magnetar/SNR assoeiation, IE 1841-045/Kes73, is located in 
figure 4 between J171405.7-381031/CTB 37B and IE 2259-1-586/CTB 109. The age of Kes73 was 
estimated by Kumar et al. (2014), as 0.75 — 2.1 kyr (table 4). Combining this with Tc = 4.7 kyr 
of IE 1841-045 (table 4), the age discrepaney of this pair beeomes Tc/tsnr = 2.7 — 8 . These two 
assoeiations do not show large overestimation faetors of Tc/tsnr as much as IE 2259-1-586/CTB 109 
association. Thus, the three magnetar/SNR associations (including IE 2259-1-586/CTB 109) suggest 
that the age over-estimation faetor, Tc/tsnr, inereases towards older objeets. This agrees, at least 
qualitatively, with the theoretieal behavior seen in figure 3 (c), as long as Pq is negligible. 

We henee tried to explain the data points of these three magnetar/SNR assoeiations with a 
eommon set of parameters, and derive a plausible range of a. Eor this purpose, three evolution tracks 
representing the solutions to equation (12) for IE 12259-1-586/CTB 109 are additionally plotted on 
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figure 4. Each parameter set is the same as that of figure 3(f). If a is small as 0 < a < 0.6 (dashed 
red line), the CXOU J171405.7-381031/CTB 37B association cannot be explained. On the other 
hand, large value of a{> 1.5) fail to explain the IE 1841-045/Kes73 association. Thus, the three 
magnetar/SNR pairs in figure 4 can be explained in a unified way if they have a common value of a 
(and also of Td) that is in the range of 0.6 < a < 1.4. 


10*^ 

10' 


lO'* 


10 ^ 


10 ^ 10 ^ 10 " 10 ^ 

'’’sNR (yr) 

Fig. 4. Relations between tsnr and Tc of single NSs associated with SNRs. Red, magenta and blue 
represent magnetars, high-B pulsars and rotation powered pulsars, respectively. Parameters are listed 
in table 4. The SGR 0501+4516/G160.9+2.6 pair is parenthesized, because the association is rather 
doubtful, and this SNR might be associated to another pulsar PSR B0458+46 (e.g., Leahy & Roger 
1991). The red, green, and blue dashed curves indicate solutions to equation (12), with (a, Td, Bq) = 
(0.6, 2.5 X 10^ yr, 6.5 x 10"^ G), (1.0, 9.2 x 10^ yr, 9.4 x lO"'^ G) and (1.4, 1.6 x 10^ yr, 
1.8 X 10^® G), respectively. They all assume Pq = 3 ms, and Bq which is specified by figure 3 (e). 
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7.2. Supporting Evidence 

The scenario so far developed implies that magnetars form a population that is much younger 
than previously thought. This important inference is supported by an independent piece of evidence. 
Eigure 5 shows a spatial distribution of NSs including magnetars. Because of steady motions after 
kick velocities are given by explosions, older pulsars with larger Tc are thus distributed up to farther 
distances from the Galactic plane. In contrast, magnetars are much more concentrated to the plane for 
their nominal age, as better seen in figure 5 (b) which is a projection of figure 5 (a) along the direction 
perpendicular to the Galactic plane. This implies two possible scenarios; magnetars, as we have 
shown, are much younger than indicated by their Tc, or their kick velocities are systematically lower 
than those of others. Recently, proper motions of four magnetars (SGR 1806—20, SGR 1900+14, 
IE 2259+586 and 4U 0142+61) were successfully measured by Tendulkar et al. (2012, 2013). They 
calculated the mean and standard deviation of their ejection velocities as 200 km and 90 km s“^. 
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respectively. They also conclude that the weighted average velocity of magnetars is in good agreement 
with the tangential velocities of the pulsar population (Hobbs et ah, 2005). Therefore, we are left with 
the former of the two possibilities. In other words, magnetars should be systematically younger than 
ordinary pulsars that have similar Tc. 


Number/Century 



D sin(6) (kpc) 


Fig. 5. (a) Spatial distributions of magnetars (red) and radio pulsars (blue). Abscissa and ordinate means distance 
from the Galactic plane and characteristic age, respectively, (b) Projection of panel (a) onto the direction per¬ 
pendicular to the Galactic plane. Radio pulsars are divided into three subgroups according to their age. (c) Age 
distributions of the objects, produced by projecting panel (a) onto the time axis. Histograms represent numbers of 
pulsars with ages in that logarithmic interval, while crosses tied by a dotted line show the object number per century. 


7. 3. Implication for the Magnetar Population 

Observationally, magnetars are no longer a minority of NS species. This is shown in figure 5 
(c), which is the projection of figure 5 (a) onto the time axis. Thus, even if Tc is not corrected for the 
over estimation, magnetars already occupy a considerable fraction of young NSs. If we replace of 
magnetars with their true ages, their dominance among young NSs will become even more enhanced. 

As yet another important implication, we expect that numerous aged magnetar descendants 
would lurk in our Galaxy. Radio pulsars, observed as the major population of NSs, cannot harbor such 
aged magnetars, because their P are shorter than those of magnetars. Instead, such objects may be be¬ 
ing discovered as weak-field SGRs, including SGR 0418-1-5729 (Rea et ah, 2013), Swift 11822.31606 
(Rea et al. 2012; Scholz et al. 2012) and 3XMM J185246.6-I-003317 (Zhou et ah, 2014). 
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8. Summary 

We performed 4 pointings observations of CTB 109 with Suzaku. The speetra extraeted from 
eastern parts of the SNR were well fitted with two plasma components having two different tem¬ 
peratures. Assuming thermal equilibrium between electrons and protons, the shock velocity was 
calculated as 460 km s“^, and the age of the SNR was estimated as 14 kyr using the Sedov-similarity 
solution. These results are consistent with the conclusion of the previous work by Sasaki et al. (2013). 
We thus reconfirmed the huge discrepancy between the age of CTB 109 and the characteristic age of 
IE 2259-1-586. 

We consider that the characteristic age of IE 2259-1-586 is significantly overestimated, as com¬ 
pared to its true age which we identify with that of CTB 109. This effect, seen also in some other 
magnetars to a lesser extent, can be attributed to decay of their magnetic fields, as implied by the 
basic concept of “magnetars”. In fact, the observed pulse period and its derivative of IE 2259-1-586 
has been explained successfully by a family of solutions to a simple equation describing the magnetic 
field decay. Eurthermore, the Tc vs tsnr relations of the three magnetar-SNR associations, including 
the IE 2259-1-586/CTB 109 pair, can be explained consistently if they have a common value of a 
in the range of 0.6-1.4. As a result, magnetars are considered to be much younger than was con¬ 
sidered so far, and are rather dominant among new-born NSs. The youth of magnetars is supported 
independently by their much stronger concentration along the Galactic plane than ordinary pulsars. 
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Table 4. Parameters for Figure 4. 


# 

Pulsar/SNR 

P (ms) 

P (ss 

B(x1012 G) 

Tc (kyr) 

TSNR (kyr) 

1 

IE 1841-045/Kes73 

11778 

4.5 X 

10-11 

730 

4 

0.75-2.1 

2 

SGR 0501+4516/G160.9+2.6 

5762 

5.8 X 

10-12 

190 

16 

4-7 

3 

J171405.7-381031/CTB37B 

3824 

5.9 X 

10-11 

480 

0.96 

0.65^253 

4 

IE 2259+586/CTB 109* 

6979 

4.8 X 

10-13 

58 

230 

10- 16 

5 

J1846-0258/Kes75 

326 

7.1 X 

10-12 

49 

0.7 

0.9-4.3 

6 

J1119-6127/G292.2-0.5 

407 

4.0 X 

10-12 

41 

1.6 

4.2-7.1 

7 

J1124-5916/G292.0+1.8 

135 

7.5 X 

10-13 

10 

2.9 

2.93 - 3.05 

8 

J1513-5908/G320.4-1.2 

151 

1.5 X 

10-12 

15 

1.6 

1.9 

9 

J0007+7303/G119.5+10.2 

315 

3.6 X 

10-13 

11 

14 

13 

10 

J1930+1852/G54.1+0.3 

136 

7.5 X 

10-13 

10 

2.9 

2.5-3.3 

11 

J1856+0113/W44 

267 

2.1 X 

10-13 

7.5 

20 

6-29 

12 

J0633+0632/G205.5+0.5 

297 

8.0 X 

10-14 

4.9 

60 

30-150 

13 

Crab 

33 

4.2 X 

10-13 

3.8 

1.2 

0.959 

14 

J0205+6449/3C 58 

65 

1.9 X 

10-13 

3.6 

5 

1-7 

15 

J1833-1034/G21.5-0.9 

61 

2.0 X 

10-13 

3.6 

5 

0.72- 1.07 

16 

J1747-2809/G0.9+0.1 

52 

1.6 X 

10-13 

2.9 

5 

1.9 

17 

J1813-1749/G12.8-0.0 

44 

1.3 X 

10-13 

2.4 

6 

1.2 

18 

J1811-1925/G11.2-0.3 

64 

4.4 X 

10-14 

1.7 

2.3 

0.96 - 3.4 


Data for P and P of pulsar were collected from ATNF Pulsar catalogue (Manchester et al., 2005)* 

Data for tsnr were collected from FeiTand & Safi-Hai‘b (2012) ^ 
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